Binding energy and fission of the heavy charged massive particle — nucleus bound state 
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We consider a possibility of capture of a heavy charged massive particle \~ by the nucleus leading 
to appearance of a bound state. A simple analytic formula allowing to calculate binding energies of 
the N\~ bound state for different nuclei is derived. If the binding energy is sufficiently large the 
particle \~ m stable inside the nucleus. The probabilities of the nuclear fissions for such states are 
calculated. It is shown that the bound states are more stable to a possible fission in comparison to 
the bare nucleus. This makes an observation of this hypothetical charged massive particle and the 
superheavy nuclei more probable. 
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I. INTRODUCTION 

The problem what the dark matter consists of is one 
of a hottest topics of the modern physics. In particular, 
a possibility of nonbaryonic dark matter is widely dis- 
cussed (see, e.g., [HQ]). < - >ne °^ possible models of 
the existence of nonbaryonic dark matter was discussed 
by Pospelov and Ritz in their recent paper [3J. They 
considered the direct and indirect detection signatures of 
weakly interacting massive particles (WIMP) x° with the 
mass m x o with a heavier, but nearly degenerate, charged 
state \ ■ The WIMP-nucleus interaction may be dom- 
inated by inelastic recombination process leading to the 
formation of N\~ bound state. 

Following this scenario and assuming that there is a 
heavy particle with the mass m x - and the charge equal to 
the electron charge e, we derive a simple analytic formula 
allowing to calculate binding energies of the Nx~ bound 
state for different nuclei. This formula can be applied 
for any Z and especially useful for heavy and superheavy 
nuclei, whose properties are less studied so far. 

In particular, we can consider the following inelastic 
scattering of \° with heavy nucleus leading to the capture 
of the particle \~ by the nucleus and the appearance of 
the bound state N\~ M : 



X + N^(N X -)+' 



(1) 



Knowing the binding energy Eq we can predict how stable 
is the bound state that occurred due to this /3 + -typc 
process. The condition of stability of the particle \~ 
inside the nucleus is \Eq\ > m x - — m x o. 

The appearance of the bound state leads, in turn, to 
that the probability of the tunneling through the fission 
barrier is changed. From general considerations we can 
expect that the N\~ bound state will be more stable to 
a possible fission in comparison with the bare nucleus. 
Our calculations show that for heavy nuclei the tunnel- 
ing probability of the bound state is reduced by a few 
orders of magnitude compared to the bare nuclei. Since 
the bound states are more stable to a possible fission in 



comparison to the bare nucleus, it opens new possibilities 
for the observations of such systems. 

The paper is organized as follows. In Sec. [IT] we derive 
the analytic formula and calculate the binding energies 
for different nuclei. In Sec. Mil we provide the expression 
for the tunneling probability and describe the Coulomb 
interaction potential for a deformed nucleus. Section HVl 
is devoted to the discussion of the obtained results and 
contains concluding remarks. The units (h = c = 1 and 
e 2 = a ss 1/137) are used throughout. 



II. BINDING ENERGY 

In this section we will find the binding energies of 
the x~ particle to a nucleus. Pospelov and Ritz [3] nu- 
merically found the binding energies of the state {N\~) 
assuming a Gaussian and steplike nuclear charge distri- 
bution for several elements and solving the Schrddinger 
equation with a given charge distribution inside the nu- 
cleus. We derive analytical formulas which reproduce the 
numerical results obtained in Ref. Q to reasonable accu- 
racy and allow us to calculate the binding energies for 
heavy and superheavy nuclei. 

Let us model a nucleus containing Z protons as a uni- 
formly charged sphere of radius R. The electric potential 
is 



V(r) = 



{Ze ( 3 _ l£ 
ii (2 2 K 
-f, r>R, 



< r < R, 



(2) 



where r is the distance from the center of the nucleus. 

Let us show that for heavy nuclei the wave function 
of the bound state N\~ will be localized well within the 
nucleus. In other words, the characteristic distances r of 
the localization of the wave function are smaller than the 
nuclear radius R even for light nuclei. For heavy nuclei 
r < R. 

As follows from Eq. ([2]) (upper line) the second term 
of V(r) is the potential of a three-dimensional harmonic 
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oscillator. The potential energy can be written as 
IZar 2 1 



U osc (r) 
where we define 



-MwV, 



2 R 3 



Za 
MR 3 



(3) 



(4) 



and M = m x -rriN /{m x - + tun) is the reduced mass of 
the particle X ~ an d the nucleus. We assume that the 
mass of the particle rn x - is comparable to or greater 
than the mass of the nucleus m^. If m x - 3> to at then 
M PS m N . 

For the harmonic oscillator, as follows from the virial 
theorem, the average potential and kinetic energies are 
equal and for the total energy of the ground state we 
obtain 



E osc = -u = Muo 2 {r 2 ). 



and, respectively, 



(r 2 



2Mui 



(5) 



(6) 



Here and in the following it is sufficient for our pur- 
poses to use the following approximations 



(7) 



is 



where ro ~ 1.2 fm, A is the nucleon number, and m 

the proton mass. 

Using Eqs. (0}, ©, and (JTJ) and denoting r av = yj (r 2 ), 
we arrive at the following expression for the ratio r av / R: 



2.7 



R A 1 / 3 Z 1 / i ' 



(8) 



In TableUwe present the results of calculation of r av /i? 
for different elements. As seen from the table, for the 
heavy and the superheavy nuclei r/R ~ 0.1. But even 
for the light nuclei we have r av < R. 

Taking into account the estimate of the characteristic 
distance r, we can neglect the wave function outside the 
nucleus and approximate the potential energy using the 
electric potential inside the nucleus only. This is the po- 
tential energy of a three-dimensional harmonic oscillator 
with the additional constant term — §4?. It is given by 



U(r) = - 



3 Za 
2~R 



+ -Muj 2 r 2 . 



(9) 



The eigen-values of the Schrodinger equation for the 
potential of a three-dimensional harmonic oscillator are 
well-known. In our case there is the additional constant 
term, so that the energies of this system are 



E 



3Za 
2~R 



+ oj (3/2 + n x + n y + n z ). 



(10) 



TABLE I: The ratio r av /R for different nuclei. 



( A N X ~) 


Z 


r av /R 


( B X ) 


5 


0.81 


( c x ) 


6 


0.75 


( 14 Nx") 


7 


0.69 


( 16 Ox") 


8 


0.64 


( 40 Ar X -) 


IS 


0.38 


( 74 Gex") 


32 


0.27 


( 132 Xe X ") 


54 


0.20 


( 202 Hg X -) 


80 


0.15 


( 232 Th X ") 


90 


0.14 


( 2S7 Fm X -) 


100 


0.13 


( 269 Ds X ") 


110 


0.13 



Let us apply this model to the case of a nucleus with 
atomic weight A. Using Eqs. (01 and ([7]) and writing the 
mass M in units of proton mass m p we obtain for the 
frequency: 



7.73 



Z nip fm 



A M 



MeV. 



(11) 



The energy of the ground state is given by 
Eq = 



-2.16 



Z fm 



A 1 / 3 r 



As we have already mentioned the characteristic dis- 
tance r is much smaller than R for heavy nuclei. For this 
reason the contribution of the second term in Eq. (fT2"|) is 
small for heavy nuclei. 

Taking into account Eq. ([7]) we arrive at the following 
approximate formulas for the frequency uj and the ground 
state energy: 



(13) 



u w 5.9 ^— MeV, 
A 




MeV. 



(14) 



In Table [TT] we list the frequencies and the binding en- 
ergies calculated for several elements using Eqs. (fT3"|) and 
(fbH) . As seen from the table there is a good agreement be- 
tween the energy values obtained in our simple approach 
and by Pospelov and Ritz Q. The heavier the atom the 
better agreement. A reason is that in our approach we 
assumed that r <C R and, as a result, neglected the part 
of the Coulomb potential for r > R. This approximation 
is good for the heavy elements while for the light ele- 
ments (like ll B) it can be used only for rough estimates. 
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TABLE II: The frequencies oj (in MeV) and the binding en- 
ergies (in MeV) of \~ f° r different nuclei. The results are 
compared (where available) with those obtained in Ref. ||. If 
\Eo\ exceeds the mass difference m x - — m x o, the particle X ~ 
is stable inside the nucleus. 



in the form 



( A N X -) 


Z 


lj (MeV) 


-E (MeV) 
This work Ref. [3] 


( n Bx~) 


5 


1.2 


2.3 


2.1 


( 12 Cx~) 


6 


1.2 


2.9 


2.7 


( 14 N X ") 


7 


1.1 


3.6 


3.2 


( 16 Ox") 


8 


1.0 


4.2 


3.7 


( 40 Arx") 


18 


0.63 


8.5 


8.0 


( 74 Ge X ") 


32 


0.45 


13.0 


12.5 


( 132 Xe X ") 


54 


0.33 


18.6 


18.4 


( 202 Hg X -) 


80 


0.26 


24.2 




( 232 Th X ") 


90 


0.24 


26.0 




( 257 Fm X ") 


100 


0.23 


28.0 




( 269 Ds X ~) 


110 


0.23 


30.3 





III. TUNNELING PROBABILITY 

In the work of Dzuba and Flambaum Q the effect of 
atomic electrons on nuclear fission was considered. It 
was shown that atomic electrons influence on the nuclear 
fission very insignificantly. The probability of the fission 
of the nuclei with Z ~ 100 is changed only at the level 
of 0.2%. 

In this paper we consider effect of a heavy charged 
particle, which forms the bound state, on the probability 
of the nuclear fission. Note that a real shape of the fis- 
sion barrier is complicated and accurate calculation re- 
quires the knowledge of the nuclear structure. At the 
same time the qualitative features in the structure of the 
fission barrier can be described by a simple parabolic bar- 
rier model [f| 0| . Since our goal is to make an estimate 
of the nuclear fission probability, the parabolic barrier 
model is sufficient for our purposes. 

The probability of the tunneling through the barrier 
can be written as [J, [6| 



1 + exp 2-7T 



\U B -E\ 



OJB 



(15) 



where Ub is the maximum of the potential energy and 
the potential barrier width ujb ~ 0.5 — 1 MeV @. 

In a case of spontaneous fission we can estimate the 
difference \Ub — E\ for typical energies -E as \Ub — E\ ~ 5 
Mev and therefore 2w\Ub — E\/uib ^> 1 ■ It is convenient 
to determine the probability of the spontaneous fission 
Po as 



Pq = exp — 2tt 



\Ub-E\ 

UJ B 



(16) 



Following Ref. Q we present the tunneling probability P 



SE 

P = P exp -2tt — 

UJ B 



(17) 



where SE = E^ nax — E^ lin is the difference between 
Coulomb energies of the particle X ~ in the points of max- 
imum and minimum of the nuclear energy U as a function 
of the deformation parameter. To find these Coulomb en- 
ergies we have to specify the form of the nucleus. It will 
be discussed in detail in the next section. 



A. Prolate ellipsoid 

Our subsequent calculations are based on the follow- 
ing experimental fact. The minimum and maximum of 
the Coulomb energy correspond to spheroidal deforma- 
tions of the nucleus. Let us consider, for example, pro- 
late spheroid with the minor semiaxis a and the major 
semiaxis c. If this spheroid was obtained as a result of 
deformation of the sphere with the radius R, then the 
condition of volume conservation reads as 



(18) 

(19) 
(20) 



a 2 c = R 3 . 

If the eccentricity e is defined by 
e 2 = 1 - a - 

2 

cr 

then for prolate ellipsoids 

a = i?(l-e 2 ) 1 / 6 , 

c = R(l-e 2 )~ 1/3 . (21) 

In Ref. 0] (see also @) it was introduced the parameter 
deformation of the sphere r] connected with the minor 
and major semiaxes of the prolate spheroid by a simple 
relation 

(22) 



c 



1 



1 



Using Eq. (fl"9"|) we obtain for the eccentricity 

2 = 3t?(6 - rj) 
(3 + ?7) 2 ' 



(23) 



B. Coulomb potential 

The Coulomb potential of the spheroid can be written 
as [|] 



U c (r,6) 



3 Za 
2 ~R 

K - 1 

~7 r ~ 



(e 2 - l)cos 2 6> + 



1 (K- 1 
2 



K\^ 2 Q\ — 



(24) 
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where the coefficient K is determined as 

2e 1-e 



(25) 



and is the angle between r and the axis z. 

It is easy to show that if the eccentricity e — > 0, i.e. 
a = c, then 



1 and 



if - 1 



and the potential Uq goes over to the potential of the 
uniformly charged sphere Uq 1 

E#(r,*) = — 



1 H 

2^2 



(26) 



C. Deformation parameter 



TABLE III: The energy difference SE (in MeV). The ratio of 
the tunneling probability P to the probability of spontaneous 
tunneling Po for the elements with Z ranging from 80 to 160 
is calculated for lob = 0.5 MeV. The notation y[x] means 
yxlO x . 



z 


A 


SE (MeV) 


P/Po 


80 


202 


0.315 


1.91-2] 


90 


232 


0.339 


1.4=1-2] 


100 


257 


0.364 


1.0[-2] 


110 


269 


0.394 


7.1 [-3] 


120 


300 


0.414 


5.5[-3] 


130 


325 


0.437 


4.1 [-3] 


140 


350 


0.459 


3.1 [-3] 


150 


375 


0.481 


2.4[-3] 


160 


400 


0.501 


1.8[-3] 



According to experimental investigations (see, 
e.g.; 0); m light actinides (like 236 U) the deformation 
77 changes from 0.2 in the minimum to about 0.8 in 
the last maximum of the fission barrier. In superheavy 
elements 77 ~ 0.4 in the last maximum. 

Using Eq. (1241) and supposing that in the minimum 
r\ = 0.2 and in the maximum r\ = 0.4 we can find SE = 
Emax ~ ^min • Taking into account that for heavy nuclei 
the wave function of the bound state is localized at the 
distances r <C R, for an estimate we can neglect the term 
- r 2 /R 2 in Eq. ([24]). Then 



U c (r,0) 



3 Za 
2 ~R 



K(l 



=-2\l/3 



(27) 



The change of the parameter 77 from 0.2 to 0.4 corre- 
sponds to change of the eccentricity e from 0.58 to 0.76. 
Then we obtain 



^<w)V 3 = 3(i r 2)1/3 m^ 

2 v ; 4e 1-e 



As 

1.494, e = 0.58 
1.475, e = 0.76 



and, respectively, 



E°(mm) r 
E c (max) : 



- 1.494 Za/R 1 
-1.475 Za/R J 



5Ek 0.019 



Za 
~R' 



(28) 



(29) 



IV. RESULTS AND DISCUSSION 

Using the derived expressions for the minimal and 
maximal Coulomb energies given by Eq. (|29p we can 
find the tunneling probabilities for the heavy and su- 
perheavy elements with Z ranging from 80 to 160. For 
calculation we again approximate the nuclear radius by 
R « 1.2 A 1 / 3 fm w 0.006 A 1 / 3 MeV -1 . For the super- 
heavy nuclei with Z > 120 whose nucleoli numbers A are 



not determined yet we put A w 2.5 Z. For an estimate 
we assume that lo b — 0.5 MeV @. 

The results of calculation of SE and P/Pq are listed 
in Table [TTT1 Our calculation confirmed that the (Nx~) 
bound state is more stable to a possible fission in com- 
parison with the bare nucleus. For the potential barrier 
width ljb — 0.5 MeV the tunneling probability of the 
bound state are 2-3 orders of magnitude smaller than 
that of the bare nucleus. The greater nuclear charge Z 
the larger difference between the probabilities P and Po- 
When Z changes from 80 to 160 the ratio P/Pq changes 
from 1.9 x 10~ 2 to 1.8 x 10~ 3 , i.e. it becomes 10 times 
smaller. It is not surprisingly because P ~ exp(— Se) and 
S E - Z/R. Since R - A 1 / 3 ~ Z 1 / 3 , S E increases as Z 2 / 3 
with increasing Z . 

To conclude, we have considered the process of capture 
of the heavy charged massive particle x~ by the nucleus 
leading to appearance of a bound state. We derived a 
simple analytic formula allowing to calculate binding en- 
ergies of the N\~ bound state for different nuclei. These 
energies can be calculated rather accurately for heavy 
and superheavy nuclei while for light elements the de- 
rived formula can be used for the estimates of the binding 
energies. 

We have calculated the tunneling probabilities for a 
number of the Nx~ bound states for the heavy and su- 
perheavy nuclei and showed that these states are more 
stable to a possible fission in comparison to the bare nu- 
cleus. Their tunneling probabilities are 2-3 orders of 
magnitude smaller than the tunneling probabilities of the 
bare nuclei. This result is important because it opens new 
perspectives to observe such bound states and get a new 
information about the hypothetical particle \~ an d the 
superheavy nuclei which were not observed so far due to 
their instability. 
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